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Abstract 

We present a model for the process of the growth of carbon nanotubes 
(CNTs) obtained by chemical vapour deposition in the presence of transition 
metal nanoparticles (Me-NPs) which act as a catalyst. We have deduced that 
the growth of a CNT occurs in the presence of two forces: (i) a viscous 
force, due to the surrounding hot gas, which opposes and slows down the 
growth of the CNT, and (ii) an extrusive force that causes the growth and 
that in the steady-state stage of the growth is completely balanced by the 
viscous force. We believe that it is the great decrease in free energy in the 
assembling reaction that occurs at the interface of the Me-NP catalyst that 
causes the extrusive force for the growth of a CNT. Moreover, the process of 
chemtsorption of a C 2 fragment, through the interaction of the C2—JT system 
with the 3d metal orbitals, has been considered as well as the coordination 
action of the Fe, Ni and Co metal surfaces. The structural properties of the 
Fe,CoandNi surfaces show that the ( 1, —1,0) planes ofFeandthe (1, 1, 1) 
planes of Co and Ni exhibit the symmetry and distances required to overlap 
with the lattice of a graphene sheet. This gives us information about the 
coordination mechanism responsible for assembling the CNTs. In fact, we 
show that it is possible to cleave an Me-NP in such a way as to match the 
correct symmetry and dimension of the armchair structure of a single-walled 
nanotube. The mechanism of C 2 addition at the edge of the growing CNT 
has also been considered in relation to the highest occupied molecular 
orbital-lowest unoccupied molecular orbital (HOMO-LUMO) symmetry. 
We demonstrate that the action of d orbitals of the metal atoms forming the 
Me-NP makes possible the thermally forbidden reaction, which involves the 
Ca-n" system. 

1. Introduction molecular nanowires. In particular, the properties of a single- 

walled nanotube (SWNT) are determined by the graphene 
Carbon nanotubes (CNTs) have stimulated a great deal of structure in which the carbon atoms are arranged to form a 
inieresi in the microelectronics industry because of their unique cylinder. This structure can be described in terms of a chiral 
conducting properties enabling the manufacture of devices vector C which considers the arrangement of the graphite 
such as field effect transistors (FETs) f 1-3], field emission hcxagons witn respect to the CNT axis. The chiral vector C is 
displays (FEDs) [4] or single electron transistors (SETs) 15, 6]. oomp^ejy defined by specifying a pair of two integer numbers 
This makes them one of the most promising materials in the (n m> identifVing ^ ^ components of C with respect to the 

fieW ^ 0,TO ^ c1 ^ n, «- ^ ^ ...... lattice parameters of the graphene sheet. In particular, SWNTs 

CNTs are hollow cylindrical shaped carbon molecules that . , . c . T , . 

^ a\^Z~^ k~ „f «c can be classified as armchair («,») nanotubes if n = m,aszig- 



because of their nanometric diameter can be thought of as 
1 Authors to whom any correspondence should be addressed. 
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zag (n, 0) nanotubes if m = 0 and as chiral nanotubes for the 
Authors to whom any correspondence should be addressed. other cases. 
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Figure 1. Weight gain versus reaction time (from [1]). Almost half 
of the growth occurs during the first minute independently on the Fe 
catalyst used, while in the following 44 min the growth occurs at a 
slower or stationary rate. 

CNTs have been obtained using several techniques such 
as arc discharge [7], laser ablation [8] and chemical vapour 
deposition (CVD) [9- LI]. In the first two cases. CNTs are 
obtained after the condensation of a hot carbon gas, while in the 
latter case the growth of a CNT is determined by the presence 
of a catalyst, usually a transition metal such as Fe, Co or Ni, 
which causes the catalytic dehydrogenation of hydrocarbons 
and consequently the formation of a CNT. In particular, the 
use of the CVD technique in the growth of CNTs is of great 
interest since it fits in with the integration process used in the 
microelectronic industry. 

These metals (Fe, Co or Ni), deposited for example on 
a nano- patterned support such as alumina or a nanoporous 
medium (porous silicon), undergo a calcination process that 
forms metal nanoparticlcs (Me-NPs). After that, the Me- 
NP catalyst is involved in a pyrolysis that, starting from a 
hydrocarbon (CH4, C2H2, CeH^X determines the growth of 
CNTs. Carbon residues are formed as a consequence of the 
catalytic dehydrogenation of the hydrocarbon molecules on the 
Mc-NP surface. 

The growth mechanisms are essentially divided into base 
growth or tip growth, depending on whether the Me-NP 
remains anchored at the support or not, respectively. This 
clearly depends on the strength of the interaction between the 
nanoparticle and the support. 

In this paper we present a model which gives clues to the 
understanding of the growth mechanism of CNTs. By focusing 
on the case of S WNTs, we will be able to consider the action 
of transition Me-NPs. And in particular, by considering their 
structural properties for the cases of Fe, Ni and Co we will 
show a correlation between armchair nanotubes and Me-NP 
structure. 

2. Extrusive-diffusive model 

A macroscopic and reasonable growth mechanism for CNTs 
can be deduced by following the development of the relative 
weight gain of the CNTs (proportional to the average length of 
the CNTs) with the reaction time as reported by Cassell et al 
[9] for the case of SWNTs and reproduced in figure 1 . From 
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Figure 2. Sketch of the growth process. The viscous force should 
be applied to the centre of mass of the system. 

the graph, we can observe that almost half of the growth occurs 
during the first minute, independently of the Fe catalyst used, 
while in the following 44 min the growth proceeds at a slower 
or stationary rate. 

Bearing this in mind, and considering that the continuous 
feedstock of carbon atoms comes from a hot and dense gas 
surrounding the growing CNT (necessary in order to gather 
many carbon atoms on the Me-NP), we can conclude that 
the CNT growth process occurs in a diffusive regime where 
surrounding hot gases provide a viscous force that slows 
down CNT growth. Moreover, in our empirical model all the 
processes that oppose the growth are also included within the 
friction force. However, one friction force alone cannot explain 
the growth mechanism. We need to invoke the presence of 
another force in order to form a CNT from an Me-NP. 

A simple model of the growth process is schematically 
depicted in figure 2, by considering that the average length 
/ of a tube is directly connected with the average growing 
speed (i.e. v = dl/dt), that the friction force can be modelled 
proportionally to the growing speed v (i.e. F — — bv where 
b is a friction coefficient) and that the extruding force F^ is 
a constant force. Therefore, according to the second law of 
dynamics: 

dv 

m— = F ex -bv (1) 

where, in the case of growth from the base, the mass m is the 
growing mass of a CNT. This equation can describe reasonably 
well the general behaviour of CNT growth in the stationary 
stage, where the stationary speed is v 3t — F ei /b. On the 
other hand, the description of the growth fails for the initial 
stage because of the discrete nature of the addition of atoms 
that cannot be described well by a continuous model. For this 
reason, near the steady state the growing mass of the CNT 
can be approximated by m = kvt + m 0 , where X is the linear 
density of the growing nanotube. Clearly, when the steady 
state is reached, the speed v is quite slow. Therefore we can 
approximate the mass of the CNT with a fictitious constant 
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fa We 1. Fitting parameters from data of HI (fitting equation 
wg = Ai + (B- A)C [1 - cxp(-//OD. 
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mass m. This allows us to integrate the preceding equation 
exactly, and by introducing the initial condition that the tube 
starts its growth with a certain speed vq wc obtain that the 
average length and hence the weight gain (wg) must follow 
the law: 

wg oc / = v s ,t -h (no - Vsi)m/b[} - exp(-6/m)f ]. (2) 

A fit of the experimental data has been tested in figure 1, 
applying this law with three parameters. Two of them arc 
proportional to v st and uo due to the wg, with the same constant 
of proportionality, while the third parameter is mfb\ the values 
of these parameters arc given in table 1. 

We observe that as uo » v Jt for all cases the growth slows 
down considerably towards the steady state. In this regime 
the model states that the viscous and the extruding forces are 
completely balanced. This empirical model allows us to follow 
the growth without any distinction between a first stage and a 
steady stage in nanotube formation and invoices the presence 
of a unique cause, i.e. the extruding force, as being responsible 
for the growth. 

What, therefore, is the origin of the extruding force that 
drives and pushes off the CNTs from the surface? In a growing 
nanotube a continuous reaction occurs at the interface of the 
Me-NP. Because of its catalytic action the Me-NP is able 
to assemble the dehydrogenated carbon fragments in a well- 
formed CNT. CNTs are stable molecules whose synthesis gives 
rise to a large and negative variation of the free energy AG. 
Indeed, we have calculated the variation of the free energy 
AG in the reaction in which a dimer is added to a molecule as 
shown in figure 3, similar to the reaction that should occur in 
the growth of a CNT. 

We estimate a variation AG = —218 kcalmol" 1 
(~10 eV/atom) at 1000 K. This huge free energy drop will be 
lost if CNTs are formed from a hot carbon gas (arc discharge, 
laser ablation), but can be used in the case of CNTs formed 
by CVD in order to do the mechanical work which is able 
to extrude the CNT. In fact AG provides an estimation of 
the maximal mechanical work that can be performed by the 
system. The Me-NP therefore favours the chemical assembly 
of carbon fragments but at the same rime directs their growth 
by determining the energy release which occurs in a well- 
defined direction. Hence, the extruding force connects with 
the — AG per unit length of the tube. Our simple model states 
the importance of the interaction at the interface between the 
carbon fragments and the Me-NP surface. 

3. Structural correlation and molecular mechanisms 

Growth models of CNTs formed by arc discharge consider 
that a CNT grows by addition of C2 (i.e. dimers) for the case 




T{K) 

Figure 3. Free energy variation for addition of C2 at the edge of the 
structure as a function of temperature. 

of armchair CNTs, by the addition of C3 (i.e. trimers) for the 
case of zig-zag CNTs, and that growth occurs at the open lip 
of the tube [12]. Similarly, in the case of SWNTs formed by 
laser ablation a spectroscopic investigation of the molecular 
species formed in the plume gives evidence for the existence 
of excited C 2 dimers [13J. 

The growth mechanism of CNTs obtained from catalytic 
nanoparticles take into account the diffusion [1 0, 1 4] of carbon 
atoms from the surface of the support into the Me-NP and the 
consequent extrusion of the walls (for multi walled nanotubes 
(MWNTs)) of the formed CNT, with the coordinating action 
of transition metals [15—18] and with a direct correlation of 
the tube diameter with that of the Me-NP [19]. Nevertheless, 
in the case of SWNTs obtained by catalytic pyrolysis of 
hydrocarbons, diffusion inside the Me-NP is not an essential 
process, whereas the interaction of carbon atoms or carbon 
fragments with the Me-NP surface plays a prominent role in the 
assembly of C atoms in an assigned CNT. In fact, hydrocarbons 
get rid of their hydrogens, eventually breaking some of their 
C bonds, and assemble on the Me-NP in order to form a CNT. 
Therefore, the chemisorption process underlies the catalytic 
action of Me-NPs. Since SWNTs have been obtained by 3d 
transition Me-NPs such as Fe, Co and Ni and since light carbon 
fragments such as C2 are involved in CNT formation, we will 
consider that a transition metal can interact with a C2 fragment 

Chemisorption of carbon fragments, such as a C2, on a 
transition metal surface is favoured due to the presence of tt 
electrons that have the right symmetry to overlap with the 3d 
electrons. In fact, by considering the molecular orbitals of a 
C2 molecule we may observe how a d\ or d a a_yj (symmetry 
e e ) orbital overlaps with a .ft?** and 7i^ py orbital respectively 
while the antibonding overlaps with a 6 xl (as also the iz^py 
overlaps with a filled d xy orbital) (t^ symmetry), in order to 
form coordinate bonds, as shown in figure 4. 

Consequently, a C2 fragment has two possible positions 
on the surface of a transition metal, i.e. on the top of the 3d 
metal atom or between two of them. Therefore, we can explain 
the coordination mechanism responsible for the assembly of 
CNTs. In fact, in order to have a cooperative efTect in the 
chemisorption mechanism, the metal atoms on a surface must 
be in the right position in order to reproduce the lattice structure 
of a graphene layer, as shown in figure 5(a). 
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Figure 4. Levels diagram for a C 2 molecule and correlations with d orbitals of a metal atom. In order to have coordinate bonds the d^ ord, r 
must be filled, as similarly the d,2_v2 or d z 2 must be unfilled. It is worth observing thai the energy spacing between the HOMO {iri py , n<i pz ) 
and the LUMO (a lpx ) of C2 is 0.26 cV (local density approximation calculation). As a consequence, in SWNTs obtained by laser ablation 
the LUMO of C 2 is quite populated. 
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Figure 5. (a) Coordination of chemisorbed C2 fragments on a metal plane reproducing the lattice structure of a graphene sheet. The thick 
black I ines represent the jr bonds of the C 2 fragments, (b) Matching between a Fe (1 , — 1 , 0), a Ni (1 , 1 , 1 ) and a Co ( 1 , 1 . I ) plane with a 
graphene sheet Fe (bec) lattice constant is a — 2.86 A, the least distance between two Fe atoms on a (1, — 1, 0) plane is d =3 2.48 A. For Co 
(1,1, l),fl = 3.54 A, d = 2.50 A. For Ni (1 % I, \\a = 3.52 A, d = 2.49 A. The lattice parameter a x a a 2 for a graphene plane amounts to 
0, = 2.45 A. 



A study of Fe, Co and Ni surfaces has shown that the ( 1 , 
—1,0) plane of Fe (bec lattice) and the (1 , 1,1) planes of Co and 
Ni (fee lattice) have the appropriate symmetry and distances to 
overlap with the lattice of a graphene sheet (figure 5(b)). This 
structural correlation suggests that if the same pattern of metal 
atoms could be transferred on an opportunely shaped Me-NP, 
where a well defined direction of growth exists, the growth of 
SWNTs of well-defined chiraiity would be possible. 

Because of this correspondence between the Fe, Co and 
Ni surfaces and the graphene layer, we have investigated the 
possibility that by cleaving an Fe crystal (and also a Co and 
Ni crystal) it is possible to shape an Me-NP in such a way that 
it matches the dimensions and symmetry of a CNT structure. 
The result of such an exercise is shown in figure 6 for the 
case of (1 0, 10) and (5, 5) armchair nanotubes for Fe and Ni 
(similar results are obtained in the case of Co Me-NPs, not 
shown here). This structural correspondence suggests hence 
that a strong correlation exists between the surface of Me-NPs 



and the SWNT structure. Clearly we suppose that the growth 
we are considering occurs at "low temperature' so that the Me- 
NPs are not in the liquid phase. 

Addition of C2 at the edge of the growing nanotubes 
is a reasonable mechanism of growth in the case of an 
armchair nanotube, because of the Ce ring closure [20]. 
This mechanism suggests that a concerted reaction can 
underlie the growth of an SWNT. For this reason the highest 
occupied molecular orbital— lowest unoccupied molecular 
orbital (HOMO-LUMO) symmetry (figure 7(a)) at the edge of 
the growing nanotube has been compared with the C2 HOMO- 
LUMO orbitals (figure 7(b)). 

In the concerted addition of the C2 to the edge of an 
armchair nanotube, during the reaction event, two elements of 
symmetry are kept: the plane perpendicular to the armchair 
structure containing the axis of the tube and the plane 
perpendicular to the axis containing the atoms at the edge. The 
HOMO-LUMO or frontier molecular orbitals (figure 7(a)) of 
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Figure 6w Me-NPs (nanorods) cleaved from the Fe (bec) and Ni 
(fee) metal structures respectively, and matching the symmetry of 
armchair (10,10) and (5, 5) nanotubes (in yellow). The metal 
surface atoms have been highlighted for clarity. The diameters of 
Me-NPs so obtained match with those of (10. 10) nanotubes 
(d = 1.4 am) and of (5, 5) nanotubes (d = 0.7 nm), rcspectiveiy. 

an 'open* nanotubc are determined by the linear combination 
of the p electrons (pj, orbitals according to the reference frame 
designated in figure 7(b)) of those carbon atoms that are at the 
edge. The symmetry properties of the two frontier molecular 
orbitals determine the reactivity of a CNT. In particular, we 
observe that the HOMO structure {figure 7(a)) is symmetrical 
(S) with respect to the plane perpendicular to the armchair 
structure, whereas the LUMO structure is antisymmetric al (A). 



In order to understand how the edge structure of a nanotube 
can interact with a C2 fragment in a concerted reaction wc 
have to consider how the HOMO-LUMO structure of the C2 
light fragments can overlap with Ihc HOMO-LUMO of the 
CNT. By comparing the molecular orbitals of C2 (reported 
in figure 4) we may observe that the <n P x LUMO of a C2 is 
symmetrical (S) with respect to the planes described above 
and can be overlapped with the HOMO of the CNT. However, 
the LUMO of the CNT does not overlap with the HOMO of 
a C2 formed by two perpendicular it orbitals (figure 4). In 
order to have a favourable overlap of the C 2 -Jr system with 
the HOMO-LUMO structure of a CNT, we need to involve an 
excited level. Nevertheless, in the case of CNTs obtained by 
laser ablation the C2 fragments are excited [13], meaning that 
the ix system must be considered in the concerted reaction. 

The role of the n orbitals of the C 2 in the 7iC 2 -CNT 
interaction can be summarized in a schematic correlation 
diagram shown in figure 7(b). In this diagram the states 
that are involved before the addition of C 2 {left) and after the 
reactive event (right) are reported and classified as symmetric 
(S) and antisymmetric (A) with respect to the two planes 
of symmetry described above. Moreover in the reaction, 
which considers that the jt bonds of a C2 fragment break and 
transform into a bonds, those states of the same symmetry 
are correlated. All the states involved in the addition of C 2 
are classified by means of the two indices S and A. The first 
index determines the symmetry of the orbitals with respect 
to the plane containing the tube axis and perpendicular to 
the armchair structure. The second index fixes the symmetry 
with respect to the plane perpendicular to the tube axis. The 
correlated levels in the absence of a catalyst are marked in 
bold. In this case the upper bonding level of the reactants 
(C 2 + edge of a CNT), which is an SA bonding state (left), 
evolves into the lowest SA antibonding state of Ihe product 
(right), determining a thermally forbidden reaction. We 
can conclude that the involvement of the n system is not 
thermally allowed and a higher level (marked with AS in the 
diagram) is required in order to obtain the addition of C2- 
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Figure 7. (a) HOMO of a (10,10) nanotube (calculated with DMol J ) [21 J. The lobes, determined by the combination of p electrons at the 
edge of the tube, are symmetrical (S) with respect to the armchair structure. The LUMO (not shown here) is similar to the HOMO but 
an ti symmetrical (A) with respect to the same plane, (b) Schematic correlation diagram of the 7rC2-CNT system for the addition of Ci to the 
armchair edge structure of the CNT involving the HOMO-LUMO structure of a (10, 1 0) nanotube and the bonding and antibonding orbitals 
formed by the Ci-n system, where the action of 3d transition metal orbitals is also considered. 
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On the other hand we have seen that the presence of a 
transition metal in the form of Me-KP acts as a catalyst and 
can help the chemisorption process and hence the coordination 
mechanism, because of the overlap of 3d orbitals with the C2 
fragments (figure 4). The same coordination of 3d orbitals 
can explain how the atoms of an Me-NP act and make possible 
thermally forbidden concerted reactions. In fact, as also shown 
in figure 7(b), an unfilled d^-^a orbital can overlap with the 
SA orbital system and lower the energy of this level (the arrow 
in figure 7(b)) describes this effect and leads to the new level) 
whereas in the same way a filled d xy orbital can overlap with 
the AS orbital system. In this way a transition from a bonding 
level system of the rcactants towards a bonding level systems 
of the products with a lower energy can occur and the reaction 
becomes thermally allowed. 

4. Conclusions 

In conclusion we have deduced that the CVD growth of CNTs 
is determined by an extrusive force opposing a viscous force. 
This extrusive force originates from the free energy release 
that occurs at the surface of an Me-NP when the C fragments 
assemble to form a CNT. This indicates the importance of the 
interaction of C fragments with the Me-NP surface. 

Chemisorpiion of C 2 molecules due to interaction with the 
d orbitals of transition metals, along with a structural matching 
of Fe, Co and Ni planes with a graphene sheet, suggest a 
mechanism for C2 formation. In particular, Me-NPs match 
the symmetry of SWNTs in the case of armchair nanotubes. A 
concerted reaction determines the assembly of the tubes, which 
is aided by the coordination action of the Me-NP catalyst. 
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